Macroscopic Coherent Rectification in Andreev Interferometers 



Jonathan Mcair 1 and Philippe Jacquod 1,2,3 
1 Physics Department, University of Arizona, Tucson, AZ 85721, 



USA 



College of Optical Sciences, University of Arizona, Tucson, AZ 85721, USA 
3 Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzerland 

We investigate nonlinear transport through quantum coherent metallic conductors contacted to 
superconducting components. We find that in certain geometries, the presence of superconductivity 
generates a large, finite-average rectification effect. Specializing to Andreev interferometers, we 
show that the direction and magnitude of rectification can be controlled by a magnetic flux tuning 
the superconducting phase difference at two contacts. In particular, this results in the breakdown 
of an Onsager reciprocity relation at finite bias. The rectification current is macroscopic in that it 
scales with the linear conductance, and we find that it exceeds 5% of the linear current at sub-gap 
biases of few tens of fieV's. 

PACS numbers: 74.45. +c, 74.78. Na, 73.23.-b 



The presence of superconductivity often magnifies 
quantum coherent effects in transport. Examples in- 
clude Aharonov-Bohm oscillations in the conductance [H- 
0] and in the thermopower coherent backscatter- 

ing [ToL 1 1 1 and resonant tunneling 12 1. The mechanism 
behind this enhancement can be traced back to Andreev 



reflection [l3j which generates new (diffuson-like) contri- 
butions to the transmission, that are sensitive to different 
phases in the superconducting order parameter or to ex- 
ternal magnetic fluxes 0, |~~ 



141. These contributions 



are proportional to the number N of transmission chan- 
nels. In purely metallic systems, quantum coherent ef- 
fects are of order one or smaller, they are therefore negli- 
gible in the limit TV ^ 1 of large linear conductances [l5j . 

Novel quantum coherent effects in transport have re- 
cently been uncovered in the form of nonlinear contribu- 
tions to the current-voltage characteristic. Of particular 
interest are contributions that are odd in a magnetic field 

b, / nl = g^{B)v 2 with gW(B) = -g&)(-B) UMI 



They originate from electronic interactions which, un- 
der finite applied biases, modify the local potential land- 
scape inside the conductor. The associated rectification 
coefficient oc dsTij has been found to be propor- 
tional to the energy derivative of a transmission coeffi- 
cient Tij, and in metallic quantum dots, it is accordingly 
sample-dependent, with a vanishing average and fluctu- 
ations decreasing with N, var(G < - 2 )) oc A^ -2 [l6, 17 [. Be- 
cause G( 2 ' is odd in B, its presence results in the break- 
down of an Onsager reciprocity relation pjj at finite bias, 
I(B,V) ^ I(-B,V). According to Mott's law, at low 
temperature the thermopower is also proportional to the 
energy derivative of the transmission |23[. Therefore, the 
question that naturally arises is whether the enhance- 
ment of the thermopower observed in mesoscopic sys- 
tems contacted to superconducting islands trans- 
lates into a similar magnification of nonlinear rectifying 
contributions to the conductance. This is the problem 
we focus on in this manuscript. 

We investigate weakly nonlinear transport in coherent 



metals connected to superconducting contacts. We find 
that the presence of superconductivity renders the rec- 
tification current finite on average and macroscopically 
large - in the sense that it scales with the linear con- 
ductance (G {2) ) = (D(N). The emergence of a finite 5 (2) 
does not require to break time-reversal symmetry in the 
metallic part of the system. It takes place, for instance, 
for two superconducting contacts with phase difference 
4> 7^ 0, 7r, {g^ 1 ) oc sin cf>. The physics behind this effect is 
that, in Andreev systems, finite biases not only modify 
the local potential landscape in the metal [24, [25 1, they 
also affect the electrochemical potential ^ sc of the super- 
conductor. In the case of a superconducting island, /x sc 
adjusts itself to ensure current conservation, and there- 
fore transmission coefficients Tij(E,eU(r), fi sc ) now de- 
pend on the absolute energy E of the charge carriers, the 
local potential landscape U(r) in the metal and addition- 
ally on ^ sc . Our key observation is that rather generic 
hybrid systems can be devised where Andreev reflection 
results in a large, finite-average derivative of with re- 
spect to the quasiparticle excitation energy e — E — /Lt sc , 
(d e Tij) = O(N) x sincf). These contributions are simi- 
lar to those giving a finite-average thermopower in An- 
dreev interferometers 0, [TH . The theory we are about 
to present predicts maximal average rectification currents 
amounting to 5-10% of the linear current at still mod- 
erate, sub-gap biases of 10-30 /iV, and for which super- 
conducting correlations persist over distances of several 
microns. In purely metallic systems, fluctuating rectifi- 
cation effects on the order of 2% typically occur for biases 
in the range 0.1-1 mV [ljH21|. 



We consider two models of Andreev interferometers, 
where two metallic terminals indexed i — L, R are con- 
nected to mesoscopic (either chaotic ballistic or disor- 
dered diffusive) quantum dots via leads carrying Ni^> 1 
transmission channels. The dots have no particular spa- 
tial symmetry and are ideally connected to two s-wave su- 
perconducting contacts with pair potentials Ae 8 ^, each 
carrying Nsi channels. Physical properties depending 
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FIG. 1: (Color online) Sketch of the two Andreev interfer- 
ometers we consider: (a) asymmetric single-cavity, and (b) 
double-cavity interferometer. The red half circles represent 
the contacts to superconductors, whose pair potentials are 
indicated. 



only on phase differences, we set 4>h = ~4>/ c 2 an d 0r = 
(j>/2 with A G 72. We consider a single superconducting 
island with two contacts into which no current flows on 
time average in steady-state. The models are sketched 
in Fig. [TJ We consider the regime where the tempera- 
ture is much smaller than the pair potential, the latter 
being in its turn much smaller than the Fermi energy, 
T <C A <C E-p- At low bias, eV <C A, the quasiparticle 
excitation energy is then always much smaller than A. 
Accordingly, we assume perfect Andreev reflection at the 
superconducting contacts. 

Both our models are specifically devised to correlate 
the average time an electron takes on its way from one 
lead to a superconducting contact, with the phase at that 
contact. This is achieved by the introduction of ballis- 
tic necks, which quasiparticles at the Fermi level cross 
in a time St. These necks are indicated by dashed lines 
in Fig. [T] The way action and Andreev reflection phases 
are correlated is easy to see by considering electrons at an 
excitation energy e, injected from the left terminal and 
Andreev reflected back to it. From Fig. [T]we see that, 
if Andreev reflection occurs at the right superconduct- 
ing contact, these electrons acquire a total phase that is 
larger by an amount 2eSr — (f> than if they hit the left su- 
perconducting contact. Such correlations were shown in 
Ref. [i| to generate large, finite-average contribution to 
the thermopower for finite cf>. We show below that they 
also result in a finite-average rectification. 

The starting point of our analysis is the scattering the- 
ory formula for the electric current in terminal i 26[ (we 
set h, fee = 1 and express electric currents in units of 
2e 2 /h) 



(i) 



with quasiparticle indices a, ft — e, +1 for electrons and 
h, — 1 for holes, the Fermi function for a /3-quasiparticle 
//(e) = (exp{[e-/3e(^ - V sc )} /T} + l)' 1 and the 
positive-defined quasiparticle excitation energy, e — \E — 
eV sc \ = \E — n sc \. We introduced transmission coeffi- 



cients for a /3-quasiparticle injected from lead j to 
an a-quasiparticlc exiting in lead i. In the presence of 
superconductivity, transmissions will depend on (i) the 
energy E of the injected quasiparticle, (ii) the local po- 
tential landscape U(r) on the quantum dot, and (iii) the 
electrochemical potential ^ sc = eV sc on the superconduc- 
tor. We take ^ sc as our reference energy and express the 
transmission probabilities in terms of two energy differ- 
ences, T^(eU(r) — fi S c,s), which describes how trans- 
mission is affected by the local potential landscape and 
the quasiparticle's excitation energy. 

At low but finite bias we expand Eq. ^ to quadratic 
order in Vj — V sc and write the current as 

* = £ g^M - v sc ) + J2 G $(^ - - v sc ) .(2) 

3 j,k 

The linear, , and quadratic, G^ k conductances are 
given by 



G 



(i) 



de (-d E f) gtj 



(3a) 



Gfl = I I de (-d £ f) (ed £ b ij S jk + 2d Vk9ij ) , (3b) 



with 



9ij 



a, P 



(4a) 



(4b) 



a,/3 



The linear conductance has been calculated earlier (sec 
e.g. Refs. 0, [T3] ) , and we focus our attention on the 
nonlinear term. The second term in the parenthesis in 
Eq. (|3b[) gives the screening contribution to the nonlin- 
ear conductance, self-consistently generated by the volt- 
age biases and the Coulomb interaction [24j. It can be 
rewritten as [H [H 



( d 9i 



\dV k 



= / dr 



S (U(r) - V sc ) 



d(U(v)-V sc ) 
dV, 



(5) 



For metallic mcsoscopic cavities the derivative of the 
transmission with respect to the local potential is random 
from one ensemble member to another, furthermore, for 
large N, it is not correlated to the local potential fluc- 
tuations [l6|. This is not modified by the presence of 
superconductivity, therefore, (dv k gij) = 0. The screen- 
ing term has no average effect on the current, and we 
neglect it from now on. The first term in Eq. (I3b| is the 
bare contribution to the nonlinear conductance [24[. In 
contrast to the purely metallic case, it depends on the 
derivative of the transmission with respect to the quasi- 
particle excitation energy e counted from the chemical 
potential of the superconductor. Below we find that the 
bare term gives a dominant, finite-average contribution 



3 



to the nonlinear conductance. We thus rewrite Eq. fl2J) 



3 



with 



G 



(2) 



p OO 

/ de (-9 £ /) 
Jo 



(6) 



(7) 



Eq. (|6|) expresses electric currents through the nor- 
mal leads as a function of the superconducting chem- 
ical potential fi sc = eV sc . In steady-state, the latter 
is self-consistently determined by current conservation 
Jl = —Ir at the normal leads. The next step is thus 
to determine ^i sc and insert its value into Eq. (J6)). One 
gets an explicitly gauge invariant expression 



i L = Q (l) v + g {2) v 2 , 

with the bias voltage V = Vl — Vr and 

= (^Mr-^MlVE^ 



kl ' 



kl 



Q {2) = J2( G LR G m ~ G RR G Ll) 



E^-^AE^) 2 



(8) 



(9a) 



(9b) 



kl 



It is easily checked that the expression for Q^> reproduces 
the result of Ref . [26| . To calculate the rectification co- 
efficient we follow the trajectory-based semiclassical 
approach of Refs. an d compute the dominant con- 

tributions to leading order in the ratio -/VsL,SR./-^ r L,R of 
the total number of transport channels in the left and 
right superconducting contacts (iVgL and N§r) and in 
the normal contacts (JVl and -/Vr). The corresponding 
diagrams are shown in Fig. O 

Eqs. ([8]-[9]) apply to both interferometers shown in 
Fig. [1] however, the coefficients G^' 2 ^ are geometry- 
dependent and we now calculate them in both cases, 
starting with the asymmetric single-cavity interferome- 
ter. For the interferometer of Fig. [TJ,, we find, to leading 
order in iV LiR , -/Vsr.sl and N S l,sr/Nl,,r, 



N l Nr 



Q {2) ) = 2esin 



N L + Nr ' 

N l NrNslN S r(Nr - N h ) 
(iV L + Afc) 4 

with a thermal damping integral 

sin(2£<5T) 



I a = I de(-d s f)d e 
Jo 



1 + (2 £ td) 2 



(10a) 
la , (10b) 

• (11) 



The value of the integral I a decreases with the temper- 
ature, T, and the dwell time, To, through the cavity. 




. (b) 
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FIG. 2: (Color online) Semiclassical diagrams that give the 
dominant phase-sensitive contributions to the rectification co- 
efficient G \ to leading order in Ns/N n [jj. (a) Contribution 
to (Tif), and (b,c,d) contributions to (T* c ). Blue (red) lines 
indicate electrons (holes) trajectories, while dashed lines in- 
dicate complex-conjugated amplitudes. Normal leads are la- 
belled i,j while superconductors are labelled Sa,S/3. Con- 
tributions to (I* h ) [(7?/)] are obtained from (a) [(b,c,d)] by 
substituting e h everywhere. 



When T < t^ 1 one has I a = 2uT St 2 csch(27rT£-r). For 
a given temperature, it is largest when St ~ 0.3/T [when 
ttTSt coth(27rT(5r) = 1]. We sec from Eq. (fTtTb) that, 
when St is nonzero, a finite average rectification current 
flows. This current is odd in <fi. Both finite <$ and St 
are required for this current to occur, because they both 
arc necessary to correlate action and Andreev phases. 
This correlation is key to obtaining a finite-average d E bij 
in Eq. Q. Eq. (|10bj) further shows that when JVj = 
0(N) » 1, i=L, R, SL, SR, and for sufficiently asym- 
metric normal terminals, |ATl — Nr\ 3> 1, the rectification 
current is macroscopically large, — O(N). 

We next present our results for the double-cavity in- 
terferometer. We find, again to leading order, 



g {1) ) = n c + 

G {2) ) = 2esin 



2iVsLiVsR 
iV SL + N S r ' 
{N Sh - N S r) 2 N S lN S rN c 
(N SL + N S r) 2 N l N r 



(12a) 
h , (12b) 



where Nq <C iVj, r is the number of channels in the neck 
connecting the two cavities and the thermal damping in- 
tegral is given this time by 



h 



de (-d e f) d £ 



Im 



exp(2ieSr) 
(l + 2ier D ) 2 



(13) 



Here, for simplicity, we took the same dwell time, td, for 
both cavities. ^From Eq. (|12b[) . we see that a macroscopic 
rectification effect also occurs in this geometry - under 
similar conditions as above, i.e. that Ni = O(N) 1 



I 



for i=L, R, C, SL and SR, and sufficiently asymmetric 
superconducting contacts, \Nsl — ATg R | 3> 1 ~ ancl that 
when T < t d \ rectification effects in the two geometries 
of Fig. [T]have the same thermal damping, — la- 

We illustrate our results in Fig. [3] for the asymmetric 
single-cavity model. We first show the current as a 
function of applied bias in Fig. [3^, for <fi = and 
4> = ir/2. We see that a rectification effect of more 
than 5% occurs at <\> = ir/2 and bias voltage of 30 /xV. 
This is rendered more evident in Fig. [3}), which shows 
the relative current asymmetry [I(V) + I(— V)]/Iu n (V), 
normalized by the linear current 7ii n = Q^V as a 
function of bias voltage. We see that at still moderate 
biases (well below the superconducting gap of Al, and 
corresponding to a coherence length vf /eV ranging 
from tens to hundreds of fim for GaAs 2DEG to 3D 
metals), the rectification effect exceeds 5%. We next 
show in Fig. [3k the rectification current as a function 
of <fi for three different voltages V — 10, 20, and 30/iV. 
In contrast to the mesoscopic rectification effects in 
metallic quantum dots which are random in an applied 
magnetic field [H, Q3, EH, we see that the presence 
of superconductivity induces a regular behavior of 
(C/( 2 )) as a function of 0, with the magnitude of the 
effect increasing with bias. Finally, the damping of the 
rectification with temperature is illustrated in Fig. 

Our approach to weakly nonlinear transport is closely 
related to the one pioneered by Biittiker and Christen 
24|. One important difference is that we here took ad- 
vantage of the presence of superconductivity to Taylor- 
expand the currents in voltages measured from the super- 
conducting potential V sc . This directly enforces gauge in- 
variance at our level of approximation, where the screen- 
ing term in Eq. (|3b[) is neglected. Current conservation is 
furthermore satisfied in our treatment by unitarity of the 
scattering matrix, JY a = A/j, and by the condition 
(sclf-consistently determining V sc ) that no current enters 
the superconducting island on time average in steady 
state. In Ref. 24 1. voltages are taken from an arbitrary 
potential as there is no superconductor. In that case 
gauge invariance is only satisfied after a self-consistent 
determination of the local potential landscape U(r) and 
of the dependence of transmission coefficients on external 
voltage biases via the latter. 

Wc have presented a theory for weakly nonlinear trans- 
port in hybrid metallic/superconducting systems and 
shown that there can be a finite average O(N) rec- 
tification for such systems. We found that, in con- 
trast to purely metallic mesoscopic systems, the pres- 
ence of superconductivity generates potentially large, 
O(N), finite-average rectification effects. The latter can 
furthermore be tuned in magnitude and direction by 
an external magnetic flux. Alternatively, we note that 
this effect leads to the breakdown of an Onsager rela- 
tion, with (still in units of 2e 2 /h) I(4>, V) - I{-(j>, V) = 
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FIG. 3: (Color online) Rectification in an asymmetric single- 
cavity interferometer with A*l = 110, Nr = 10, Nsh = iVsR = 
30, T = lOOmK, Sr = 0.3/T, T < r D 1 . (a) Electric current at 
±V vs. the bias voltage \V\. At <f> = and vr, /(V) = -I(-V) 
and there is no rectification. At <j> = it/2, rectification is 
maximal (red curves). A subdominant oscillating contribu- 
tion to the linear current (see Eq.(3a) in Ref. [y|), not in- 
cluded in Eq. (|10a[) . is taken into account here, (b) Relative 
current asymmetry /(V) + I( — V) normalized with the linear 
current 7ii n (V) = G"'V for four different superconducting 
phase differences, (c) Sample-averaged rectification conduc- 
tance {Q ( - 2 ' I )V vs. <j) for three different bias voltages. The 
oscillating part of the linear conductance is shown for com- 
parison (dashed line), (d) Thermal integral I a of Eq. (|lip 
giving the damping of (G^) with temperature for the single- 
cavity model of Fig. QJi. 



4e sin^[JV L JV R JV SL JV SR (JV R - A^l)/(^l + A^r) 4 ]/^ 2 for 
the asymmetric single-cavity model and /(</>, V) — 
I(-4>,V) = 4esin0[(AT SL - A^sr) 2 A^ sl A^ sr A^c/(A^ sl + 
A r s R ) 2 N^Nn]!^ 2 for the double-cavity model. We ex- 
pect the rectification effect we predict to be experimen- 
tally testable in Andreev interferometers such as those of 
Refs. @,|,@. 
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